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Abstract 

New lattice model for the gradient elasticity is suggested. This lattice model gives 
a microstructural basis for second-order strain-gradient elasticity of continuum that is 
described by the linear elastic constitutive relation with the negative sign in front of the 
gradient. Moreover the suggested lattice model allows us to have a unified description of 
gradient models with positive and negative signs of the strain gradient terms. Possible 
generalizations of this model for the high-order gradient elasticity and three-dimensional 
case are also suggested. 


1 Introduction 

The two most widely used theories of elastic deformation in solid materials are a microscopic 
approach based on the statistical mechanics of lattices PEIS] and the quantum theory of 
solid-states pf], and a macroscopic approach based on the classical continuum mechanics PE]. 
Continuum elasticity is a phenomenological theory representing continuum limit of lattice dy¬ 
namics, where the length-scales are much larger than inter-atomic distances. Nonlocal elasticity 
theory is based on the assumption that the forces between material points can be at long-range 
in character, thus reflecting the long-range character of interatomic and intermolecular forces. 
In general, the nonlocal continuum models describe materials whose behavior at any point de¬ 
pends on the states of all other points in the media, in addition to its own state and the state 
of external fields. Such considerations are well-known in solid-state physics, where the nonlocal 
interactions between the atoms and molecules are prevalent in determining the properties of 
the media and materials. 

The theory of nonlocal continuum mechanics was formally initiated by the papers of Kroner 
P and Eringen, Edelen (U [9]. Kroner [7J indicated the relation between nonlocal elasticity 
theory of materials with long range cohesive forces. Eringen and Edelen [9] provided derivation 
of the constitutive equations for the nonlocal elasticity. Eringen and Kim [lOj described a 
relation between non-local elasticity and lattice dynamics. Kunin described the physical aspects 
of nonlocal elasticity in the book PP, and studied various problems in Fourier space. In the 
book p2] Eringen considered a unified approach to field theories for elastic solids, viscous fluids, 
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and heat-conducting electromagnetic solids and fluids that include nonlocal effects. Rogula Pi 
considered the mathematical aspects of nonlocal elasticity models, proposed different types of 
nonlocal constitutive relations between stress and strain, and applied it to various problems 
in continuum mechanics. Non-local continuum mechanics has been treated with two different 
approaches [13, OH]: the gradient elasticity theory (weak non-locality) and the integral non¬ 
local theory (strong non-locality). In this paper we discuss the gradient models of non-locality 
elasticity. Usually two classes of gradient models are distinguished by the different signs of the 
strain gradient terms in the constitutive relations for the strain and the stress a if. 
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where A and /i are the Lame coefficients, l is the scale parameter. If l 2 = 0, we have the classical 
case of the linear elastic constitutive relations for isotropic case that is the well-known Hooke’s 

law. 

The first class of gradient elasticity models are described by equations (P with the positive 
sign in front of gradient. The main motivation to use this form of the gradient elasticity is 
the description of dispersive wave propagation through heterogeneous media. In many studies, 
gradient elasticity models with the positive sign in ([Tj) have been derived from associated lattice 
models by the continualization procedure for the response of a lattice nsnnnm 

The second class of gradient elasticity models are described by equations © with the 
negative sign in front of the gradient. The strain gradients in equation (P) with the negative 
sign are equivalent to those derived from the positive-definite deformation energy density, and 
therefore these models of the strain gradients are stable. 

The positive sign of the strain gradient term in equation (JT|) makes this term destabilizing. 
The corresponding equation for the displacements is unstable for wave numbers k > l/l 2 
mmmm- In dynamics the instabilities lead to an unbounded growth of the response in 
time without external work. It is known the instabilities are related to loss of uniqueness in 
static boundary value problems. Instabilities in statics and dynamics for the second-gradient 
models with the positive sign are discussed in [21]. 

At this moment there is the opinion that gradient elasticity models with the negative sign in 
equations (P cannot be obtained from lattice models pH] • It is usually assumed that this class 
of the second-gradient models does not have a direct relationship with discrete microstructure 
and lattice models [15]. It was proved that the homogenization (continualization) procedure, 
which is considered in [P PI HH1 3,22], uniquely leads to a second-order strain gradient term 
that is preceded by a positive sign. The second-gradient model with negative sign cannot be 
derived by this homogenization procedure. From a mathematical point of view it is caused by 
properties of the Taylor series that is used in this procedure. 

In this article we propose lattice models, that allow us to derive linear elastic constitutive 
relations with negative and positive signs. Moreover the suggested lattice models give unified 
description of the gradient models with positive and negative signs of the strain gradient terms. 
To obtain continuum equation from the lattice equations we use an approach that is suggested 
in [23123 m w\. 
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2 Equations of lattice model 


Let us consider the vibration of an unbounded homogeneous lattice, such that all particles are 
displaced from its equilibrium position in one direction, and the displacement of particle is 
described by a scalar held. We consider one-dimensional lattice system of interacting particles, 
where the equation of motion of nth particle is 

d 2 u ft) + °° + °° 

m^ = 92 £ K 2 (u<, tyi) u rn {t^) ^ ^ K 4 (n,m) “I - (2) 

m=—oo m =—00 

m^n m^n 

where u n {t ) are displacements from the equilibrium, #2 and g± are coupling constants, F{u n ) 
is the external on-site force, K 2 (n,m) and K^n^m) are the functions with different power-law 
asymptotic behavior of the functions 

00 

K s (k) = 2 K s (n, 0) cos(fcn), (s = 2; 4) (3) 

n =1 

for k —> 0. We will consider interactions terms for which the difference K s {k ) — iL s (0) are 
asymptotically equivalent to |fc| s as \k\ —> 0. Note some general properties of K s (n,m), with 
s = 2; 4. The conservation law of the total momentum in the lattice ([2]) in case of absence of 
external forces Ffu n ) = 0 gives 

+oo +oo 

£ K 2 (n,m) = 0, ^ K 4 (n,m) = 0, (4) 

m =—oo m =—oo 

m^n m^n 

for all n. For homogeneous unbounded lattice, we have 


AT 2 (n, m) = A' 2 (n — m), A' 4 (n, m ) = K 4 (n — m), 

where elements of K s (n,m) are constrained by the conditions ((4|), and 

K s (n — m) — ^ K s [n — m) = 0. (5) 

m n 

For a simple case each particle is an inversion center and K s fn — m) = K s (\n — m\), where 
s = 2; 4. Using the condition 03]), we can rewrite equation 02]) as 

M—^ = g 2 K 2 (\n - m\) (u n - u m ^j + g 4 ^ K 4 (\n - m\) (u n - u m ^j + F(u n ). (6) 

m =—oo m=—oo 

m^n m^n 

In this form of equation of motion the interaction terms are translation invariant. It should be 
noted that the noninvariant terms lead to the divergences in the continuum models [2?]. 
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Let us give an effective discrete mass-spring system for the suggested lattice model ([6]). 
In Figure 1, we present the nearest-neighbor and next-nearest-neighbor interactions only. In 
general, the functions K 2 (\n—m\) and K±(\n—m\) describe long-range interactions with power- 
law asymptotic of ([3]). 


l,eff l,eff U e f f 

ft' 4 ft '4 ^4 



Figure 1: Discrete mass-spring system with effective stiffness coefficients = ^2 (< 72 , 5b) 
and k e /f = k e J*(g 2 , < 74 ), the mass M and the distance a that correspond to the lattice model 
with coupling constants g 2 and g±. 

3 From lattice model to continuum model 

Let us consider a set of operations pi, 25j [2?] that transforms the equations of motion of the 
lattice model into a continuum equation for the displacement field u(x, t ). We assume that u n (t) 
are Fourier coefficients of the field u(k,t ) on [—K 0 /2, K 0 /2] that is described by the equations 

+OO 

u(k,t) = ^2 u n (t) e~ lkXn = J 7 &{u n {t)}, (7) 

n=—oo 
| r+Ko/2 

U n (t ) = — / dk u(k,t) e lkXn = (8) 

A 0 J-Ko/2 

where x n = na and a = 2n /Kq is distance between equilibrium positions of the lattice particles. 
Equations (J7J) and (jHJ) are the basis for the Fourier transform J 7 ^ and the inverse Fourier series 
transform J 7 ^ 1 . The Fourier transform can be derived from (O and (jHJ) in the limit as a —» 0 
(Kq —» 00 ). In this limit (a —* 0 or K 0 —> 00 ) the sum becomes the integral, and equations (J7J) 
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and (JED become 


u(k,t ) = 


u(x, t ) = 


r +oo 

' —oo 
/»+oo 


2tt 


dx e lkx u(x,t) = .F{w(:e, £)}, 
dfc e ikx u{k,t) — J r ~ 1 {u(k,t)}. 


(9) 

( 10 ) 


Here we use the lattice function 


Un{t) 


2vr 


w(x n ,t) 


with continuous u(x,t), where 


27 in 


x n = n a = —-> a;. 

Ao 

We assume that u(k,t) = Cu(k,t), where £ denotes the passage to the limit a —* 0 (A'o —> oo), 
i.e. the function u(k,t ) can be derived from u(k,t) in the limit a —> 0. Note that u(k,t ) is a 
Fourier transform of the field u{x, t). The function u(k, t ) is a Fourier series transform of u n {t), 
where we can use u n {t) = (2 tt/K 0 )u(n a, t). 

We can state that a lattice model transforms into a continuum model by the combination 
T l C T a of the following operation. 

The Fourier series transform: 


A a : 

The passage to the limit a —> 0: 

£ : 

The inverse Fourier transform: 

T- 1 : 


U n {t) Aa {U n (t)} = U(k,t). 

(ii) 

u(k,t ) —> £{u(k,t)} = u{k,t). 

(12) 

u(k,t ) —> J r ~ 1 {u(k,t)} = u(x,t). 

(13) 


These operations allow us to get a continuum model from the lattice model [2T1 25 ) \27] . 


4 Lattice model with nearest-neighbor interaction 

Let us derive the usual elastic equation from the lattice model with the nearest-neighbor in¬ 
teraction with coupling constant g 2 = K by the method suggested in [231E25, |27]. We will use 
equations ([2]) with 

+oo 

E K 2(n, m ) u m (t) = u n+ i(t) - 2 u n (t) + w n -i(t), K 4 (n, m) = 0, (14) 

m =—oo 
m^n 

where the term K2(n,m) describes the nearest-neighbor interaction. 
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We can give the following statement regarding the lattice model with the nearest-neighbor 
interaction and the corresponding continuum equation that is obtained in the limit a —» 0. 
Proposition 1 . In the continuous limit (a0) the lattice equations of motion 


d 2 u n {t) 

dt 2 


K ■ 


2 U n (t) + -Un-l(t)) + F(u n (t )) 


(15) 


are transformed by the combination T l C Fa of the operations 0711731) into the continuum 
equation: 


d 2 u(x,t) _ 2 d 2 u{x,t) 1 ,, x 

dt 2 e dx 2 + 


(16) 


where C 2 = E/p = K a 2 /M is a finite parameter, and f(u) = F(u)/(Aa). 

Proof. To derive the equation for the field u(k, t ), we multiply equation (113]) by exp (—ik n a), 
and summing over n from — oo to +oo. Then 


+oo ,2 4-oo +o° 

M Y e ~ ikna ^r = K - E e- ikna (u n+1 -2u n + u n _^+ e~ ikna F(u n ). (17) 

n =—oo n=—oo n =—oo 


The first term on the right-hand side of (TT7T) is 


+oo 


+oo 


K ■ Y e~ ikna K 2 {n,m)u n = K ■ Y 


—ik n a 


^n+1 2 U n + r H n —\ ) 


+oo 


+oo 


4-oo 


K ■ Y e~ iknd u n+1 -2K ■ Y 


e~ iknd u, 


+ K- £ 


e - ikna u n _ 1 = 


n =—oo 
4-oo 


77.— — OO 

4-oo 


4-oo 


= e 


e~ ikna u, 


+ e~ ikd K ■ Y 


F ikja u, 


3 =~°o 


f ka K ■ Y e~ ikma u m -2K ■ Y 

171 = — OO 77 .— — o 

Using the definition of u(k,t), we obtain 

4 -oo 

K ■ E e~ ikna Iu(n, m)u n = K • ( e ikd u(k, t ) - 2u(k, t ) + e~ ika u(k, t )) = 


= K ■ [e lka + e lka — 2 S ju(k, t ) = 2 K ■ ^cos (ka) — lj u(k , t) = —4K ■ sin 2 (^~^j 7(U t). (18) 
Substitution of (1T8]) into (1171) gives 

m~^t^ = ~^ K ' sin2 *) + Fa{F (w n (f))}. (19) 
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For a —* 0, the asymptotic behavior of the sine is sin(k a/2) = ka /2 + 0((ka) 3 ), then 


-4 sin 2 



(k a) 2 + 0((k a) 4 ). 


Using the finite parameter U 2 = K a 2 /M, the transition to the limit a —> 0 in equation ([T9]) 
gives 

= -C 2 e k 2 u(k,t) + jjT{F(u)}, 

where 

M „ Ka E Ka 2 


( 20 ) 


P = 


A a' 


E = 


The inverse Fourier transform T 1 of 


A ’ 
has the form 


= _= 

e p M 


( 21 ) 


dt 2 


~C 2 T 1 {fc 2 h(fc,f)} + 


where /(u) = F(u)/(Aa) is the force density. Then we can use E^ l {u{kA)} = u(x,t), and the 
connection between derivatives and its Fourier transform: J r ^ 1 {k 2 u(k,t)} = d 2 u(x,t)/dx 2 . As 
a result, we obtain the continuum equation (lT6]h This ends the proof. 


As a result, we prove that the lattice equations (1T5]) in the limit a —» 0 give the continuum 
equation with derivatives of second order only. This conclusion agrees with the results of [23], 
where the relation 


exp i 



u(x,t ) 


u(x + a, t ) 


and the representation of (fT5j) by pseudo-differential equation are used. 


5 From general lattice model to gradient elasticity model 

Let us consider the lattice model that is described by (|2|) . where the terms K s (n,m ) with s = 2 
and s = 4 satisfy the conditions 

OO 

K s (n,m) = K s {\n - m\), ^ \K s (n)\ 2 < oo. (22) 

n =1 

To describe gradient elasticity models, we consider the inter-particle interactions, that are 
described by K s (n ) (s = 2 or s = 4) of the following special type. We assume that the function 

+oo OO 

K s {k)= e ~ iknR s(n) = 2j2 K s(n)cos(kn), (23) 

n=— oo n= 1 

n/0 
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satisfies the condition 


lim 

fc-> o 


Ks(k)-k a ( 0) 

\k\ s 


= A a , 


where 0 < |A S | < oo. Condition (1241) means that 

K a (k) - K a (0) = A a \k\ 9 + R a (k), 


(24) 


(25) 


for k —> 0, where lim fc ^ 0 Rs{F)/\k\ s = 0. This also means that, we can consider arbitrary 
functions K s (\n — m\) for which K s (k ) — K s { 0 ) are asymptotically equivalent to |/c| s as \k\ —> 0. 

As an example of the interaction terms K s (|n — ?n|), which give the continuum equations 
of gradient elasticity models, we consider the function 


K s (\n - m\) = 


(- 1) 1 


n—m\ 


2T(s/2 + 1 + \n — m\) T(s /2 + 1 — \n — m \)' 


(26) 


We use 2 in the denominator to cancel with 2 from equations f[23|) . The terms K s (|n — m\) 
are considered for n ^ m, i.e. \n — m\ ^ 0. For s = 2 j, we have K s (\n — m |) = 0 for all 
\n — m\ > j + 1. The function K s {n — m) with even value of s = 2 j describes an interaction 
of the n-particle with 2 j particles with numbers nil. . . n ± j. To represent properties of 
), we can consider the function 


fi<{x,y) = Re[iF 2/ (x)] = 


Re[(—1)1^1] 


2T(y/2 + l + \x\)T(y/2 + l-\x\) 


(27) 


of two continuous variables x and y > 0. Note that Re[(—1) 1^1] = (—1)^ for integer x = n — m. 
The plots of the function (|27l) are presented by Figures 2 and 3 for different ranges of x and y. 
This function decays rapidly with growth x and y. The function (1271) defines the interaction 
terms K s (\n — m |) by the equation K s (\n — m\) = //f(|n — m\, s). 

Using an inverse relation to (123|) with K s (k ) = \k\ s that has the form 


K s (ri) — — I k s cos (nk)dk 


7r 


we get another example of K a (\n — m\) in the form 


K s (\n - m\) = 


7T 


S + 1 


\F'i 


s + 1 1 s + 3 
2 ’ 2 ’ 2 ’ 


7T 2 (n — mf 


(28) 


where \F 2 is the Gauss hypergeometric function (see Chapter II in J3U]). Note that the inter¬ 
actions with (l28]i for s = 2 and s = 4 are long-range interactions of n-particle with all other 
particles (m e N). It is easy to see that expression (1281) is more complicated than fl26j) . 

For s = 2, we can also use the long-range interactions in the following two forms 











A main advantage of the interaction in the forms (1261) and ()28l) is a possibility to use for 
other generalizations for the case of the high-order gradient elasticity by using arbitrary integer 
values of s and the fractional generalization of gradient elasticity by non-integer values of s. 
Proposition 2. The lattice equations 


M 


d 2 u n 
dt 2 


+oo +oo 

g 2 y K 2 {\n - m\) [u n - u m ] + g 4 y Kfi\n - m\) [u n - u m \ + F(u n ), (30) 

m=—oo m =—oo 

m^n m^n 


where g 2 and g 4 are coupling constants, K 2 (\n — m \) and K 4 (\n — m \) are defined by TWil . 
u n = u n {t), are transformed by the combination F l C J -a of the operations 0711731) into the 
continuum equation: 


d 2 u(x,t ) d 2 u(x,t) d 4 u(x,t) 

dt 2 2 dx 2 + 4 dx 4 


1 

P 


f(u(x,t)) = 0, 


(31) 


where 

n a 2 n g±a A 
2 AM' 4 48M 

are finite parameters, p = M/(Aa ) is the mass density and f(u ) = F(u)/(Aa ) is the /orce 
density. 



Proof. To derive the equation for the field w(fc, t), we multiply equation (l30|i by exp (—ik n a), 
and summing over n from — oo to +oo. Then 


+oo 


M Y, e ~ ikUa %. u ^= Y Y Y e ^ kna 9 s K s(\n-m\) [u n -u m }+ y e~ ikna F{u n ). 

DO 

(33) 

(34) 


+oo +oo 


+oo 


n =—oo m =—oo s= 2;4 
m^n 


The left-hand side of (1331) gives 


- iknad 2u n(t ) _ ^ -ikna /,>, _ 9 2 u(k,t) 

2^ e fp2 ~ fp. 2 2^ e ~ 


dt 2 ’ 


where u[k,t) is defined by (J7|). The second term of the right-hand side of equation (l33l) is 
^2n=-oo e~ lkna F{u n ) = Fa{F(u u )}. The first term on the right-hand side of (1331) is 

+oo +oo 


y, y e 2fcria Ar s (|n - m\) [u n - u m ] = 


n =—oo m =—oo 
m^n 


+oo +oo 


—)—C50 +00 


y y e -^-A 3 (|n-m|) Mn - y y e- ikna K s {\n-m\)u m . ( 35 ) 


n=—oo m =—oo 

rriy^n 


n =—oo m=—o o 
m^n 
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Using (JTJ) and (|22lh the first term in r.h.s. of (1351) gives 


+OO +00 +OO 

—ik n a 


+oo 


53 I] e- ikna K s (\n-m\)u n = E K s(m!) = u{k,t)K s ( 0), (36) 


n=—oo m =—oo 
m^n 


m =—oo 
m'/0 


where 


-(-OO 


AVfcaj - £ e-“”“JC(n) = J- A {A',(n)}. 


(37) 


n=—oo 

n/0 


The second term in r.h.s. of (1351) gives 

+oo +oo 


+oo 

E 53 - m|) J3 


+oo 


Um. 


n =—oo m =—oo 
m^n 


n =—oo 
n^m 


m =—oo 


+oo +oo 

= 53 e~ ikn ' a K s {n') 53 u m e~ ikma = K s (ka) u(k,t). 

n'=—oo m=—o o 

nVO 

As a result, equation (1331) has the form 


m 5VM) = ^(0) _ AT s (fca)) t*(M) + ^{%)}, 

s=2;4 


(38) 


(39) 


where J r A{A(n n )} is an operator notation for the Fourier series transform of F{u n ). 
Using the series (see Sec.5.4.8.12 in [28]) of the form 


E 


-i)’ 


—' T(z/ + 1 + n)r(z/ + 1 — n) 


cos (nk) 


2 2v - 1 
T(2is + 1) 


Sill 


2u 



1 

2 r 2 (z/ + 1 )’ 


where zz > —1/2 and 0 < k < 2i r, we get for the function (f23l) of the form (1261) the equation 


K s (ak)-K s ( 0) 


2 s ” 1 

m+i) 


sin 3 



1 

2r(s +1) 


ak\ s + 0(k s+2 ). 


(40) 


Here we use v = s/2 and sin(/c/2) = k/2 + 0(k 3 ). Note that 2 in the denominator of (1261) 
cancels with 2 from equation (T23T) in front of the sum from zero to infinity. The limit k —> 0 
gives 


K s (k)-K s ( 0) 1 

fc™ \k\ s 2T(s +1) ’ 


(41) 


and we have A s = l/(2T(s + 1)). 
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The Fourier series transform J-\ of (1301) gives (1391) . We will be interested in the limit a —> 0. 
Using (l40lh equation (f39|) can be written as 

^ %Ak) U(k,t) - g -F- t 1A (k) u(k,t) - TjF a {F („„(())} = 0, (42) 

where 

UW - - 2r(s 1 + 1) |fc| J + a 2 0(|fc|«+ 2 )- (43) 

In the limit a —> 0, using 

%(k) = CT, A (k) = - * + (s = 2;4), (44) 

we get 

t I (k)=Ct IA (k) = -\\k\\ t i (k)=Ct lA (k) = -^\k\ i - (45) 

The passage to the limit a —> 0 for the third term of gives FaF( u n ) —> CJ r / \F(u n ). Then 

CF A {F(u n )} = F{CF(u n )} = F{F(Cu n )} = F{F(u(x, t))}, (46) 

where we use CFa = FC. 

As a result, equation (|42|) in the limit a —> 0 gives 

^u(k, t ) - G 2 t 2 {k) u(k , t ) - G a %(k) u(k, t ) - y -F{F (u(x, t))} = 0, (47) 

ot z M 

where u(k,t ) = Cii(k,t), and we use finite parameters G 2 and G A , that are defined by (l32]h 
The inverse Fourier transform of (|47|) is 

0 2 1 

wry u(x, t ) - G 2 T 2 (x) u(x, t) - G 4 %{x) u{x , t) -/ (u(x, t)) = 0, (48) 

ot z p 

where the finite parameters G 2 and G A are defined by (l32lh Using (148|) the operators F 2 (x) and 
7l(x) are defined by 

T 2 (x)=T- 1 {t 2 (k)}=+F T i (x) = F-\t i (k)} = -F (49) 

ffere, we have used the connection between the derivatives of the second and fourth orders and 
their Fourier transforms k 2 <—> —d 2 /dx 2 and k 4 <;—> Fd^/dx 4 , and equations (1451) . 

As a result, we obtain continuum equations (1TB . This ends the proof. 

Proposition 2 illustrates the close relation between the discrete microstructure and the 
gradient non-local continuum. Let us consider special cases of the suggested model. 


11 










Lattice equations (l30li have two parameters g 2 and g±. The corresponding equation (13B for 
the elastic continuum has two finite parameters G 2 and G 4 . If we use g 2 = 4 K and (74 = 0, 
then G 2 = Cl = K a 2 /M , G 4 = 0, and we get equation (fT 6 ]h If we assume that r / 2 = 4 K and 
<74 = —4/i, then G 2 = C 2 = A' a 2 /M, G 4 , = C 2 a 2 /48 and we get the equation 


d 2 u(x,t) _ 2 d 2 u(x,t) a 2 C 2 d 4 u(x,t) 1 
9t 2 e 9a ; 2 12 9a ; 4 p U 


(50) 


where C e = \JE / p is the elastic bar velocity. Equation (1501) can also be derived by the 
homogenization procedure mmm- 

In general, the coupling constants g 2 and g 4 are independent. Therefore the coupling con¬ 
stant g 4 may differ from the constant g 2 = 4 A'. If the relation of stress and displacement of 
the form e(x,t ) = du(x,t)/dx is used, and the continuum equation (13B is expresses as 

d 2 u(x,t ) da(x, t) 

where p = M/(Aa), then the constitutive relation can be represented by 


a = E 


g 2 _ g 4 a 2 9 2 e \ 

4 K ~ 48A' dx 2 )' 


(51) 


where we use E = K a/A. Therefore, using the correspondence principle, we will assume 
g 2 = 4 K. The second-gradient parameter l is defined by the relation 

,2 \9a\ a 2 

48 K ’ 

where the sign in front of the factor l 2 in the constitutive relation is determined by the sign of 
the coupling constant g 4 . If the constant g 4 is positive then we get the second-gradient model 
with negative sign. As a result the second-gradient model with positive and negative signs 

a = E^e- sgn(g 4 ) l 2 (53) 

can be derived from a microstructure of lattice particles by suggested approach. The proposed 
model as shown above uniquely leads to second-order strain gradient terms that are preceded by 
the positive and negative signs. It should be noted that positive value of coupling constant g 4 of 
lattice model can lead to effective stiffness coefficient of the next-nearest-neighbor interaction 
with non-convex elastic energy potentials in the effective discrete mass-spring system. The 
strain gradients in continuum equation with the negative sign are equivalent to those derived 
from the positive-definite deformation energy density, and therefore these continuum models 
are stable. The lattice models with negative value of coupling constant g^ of lattice model leads 
to the continuum equation with the positive sign in front of the parameter l 2 . This continuum 
equation is unstable for wave numbers k > l/l 2 The instability leads to an unbounded growth 
of the response in time without external work. 
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6 Possible extensions of general lattice model 

The suggested lattice model can be generalized and extended for the high-order gradient elas¬ 
ticity and for three-dimensional lattice models. Let us give some details about these general¬ 
izations. 

We can consider a generalization of the suggested lattice model by using the sum of the 
functions (1261) with the even value s. Using the functions (|26|) with s = 6 and other even values, 
we can consider the lattice models for high-order gradient elasticity [I5j 29, [22]. We can state 
that the lattice equations 

M g2 i S K v(\ n ~ m \) (u n {t) - Um(*)) + F(u n ), (54) 

j=l m= — CO 

m^n 

where g 2 j (j = 1 are coupling constants, and K 2j (|n — m\) are defined by f[26|) . are 

transformed by the combination C J -a of the operations (I111I13P into the continuum equation 


d 2 u(x,t) ^ ■ d 2j u(x,t ) 1 

dt 2 + ZJ -1 ) - ~J( U ) = °> 


3 = 1 


dx 2 i 


P 


(55) 


where G 2 j = g 2j a 2] /(2 F(2j + 1)), (j = 1 are finite parameters. The proof of this 

statement is similar to proof of Proposition 2. 

The suggested one-dimensional lattice model for second-gradient elasticity can also be gen¬ 
eralized for the three-dimensional case. We may consider a three-dimensional lattice that is 
described by the equations 

J2J 

-^T= 5^ K ?( n - m) (u l Q - u l m ) + ^f(n-m ) (u l a - u l m ) + F k (u n ), (56) 

m:m^n m:m^n 


where n = (ni,n 2 , 77 . 3 ), k, l = 1, 2, 3 and we assume a sum over repeated index l — 1, 2, 3. In the 
model (1561) the coupling constants are included in the tensors K^ l (n — m) = Kg l (m — n) that 
are distinguished by different power-law asymptotic behavior. We assume that the functions 
Kg l { k) — Kg l ( 0) are asymptotically equivalent to ki kj and kikj |k| 2 for s = 2 and s = 4 
respectively, where 

U'(k) = T e_ik “ K s‘( n )- 

n 

To get continuum equation, we consider the field u n (t) as Fourier coefficients of the function 
w(k, t), where k = (ki, k 2 , k 3 ), by 


u\k,t) = J2 <(*) e" ?:krn , 

n 
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where r(n) = r n = y^, 3 _ 1 rij a,;, with the translational vectors aj of the lattice. In three- 
dimensional lattice model for second-gradient elasticity, we should consider interaction terms 
K kl {n — m) that satisfy the conditions 


lim 

ki ,kj —>-0 


Af (k) - K kl (0) 


kj k-i 


AS ( 2 ), 


lim 

ki,kj,\ k|—>0 


Kf (k) - Kf (0) 

U. U. lb-12 
n>2 | 1 



where A k j(s) are the coupling constants. In the continuous limit (|a,; 
dimensional lattice (156]) gives the continuum equations in the form 


( i,j = 1,2,3), 

(57) 

—* 0), the three- 


d 2 u k (r , t ) 

W 2 


,<9V(r,£) kl d 4 u l (r,t) 

Gij (2) —— -h G i • (4) 7-—-— --— 

J OXiOXj J OXipXjOXm OX m 


-f k (u(r,t)) = 0, 
P 


(58) 


where we assume a sum over repeated indices i, j, l, m G {1, 2, 3}, and 




3 l II II 

nki(A\ — V' l a?; l l a jl l a?n 

~ m 

m= 1 



(59) 


where no summation over repeated indices. We can consider the case with G k J( 4) = l 2 G k j( 2), 
where G k j( 2) = Cikji can be considered as a stiffness tensor and l 2 is the scale parameter. For 
isotropic case, we have Cijki = \8ijSki + p(&ikbji + &ubjk)- In general equation (1581) describes 
anisotropic gradient continuum. A more detailed description of the three-dimensional lattice 
model (156|) will be made in the following article. 


7 Conclusion 

In this paper lattice models for strain-gradient elasticity of continuum are suggested. The 
first advantage of the suggested lattice models is a possibility to consider these models as a 
microstructural basis of unified description of gradient models with positive and negative signs 
of the strain gradient terms. A second advantage of the proposed model is that it can be easily 
generalized to the case of the high-order gradient elasticity by using the other even values of s. 
Using (126|) with positive integer s = 2j, we have non-local interaction of the lattice particle, that 
gives the derivatives of integer order 2 j in the continuum equation. Three-dimensional lattice 
models and the correspondent continuum equation can also be formulated as ()56l) and ()58lh 
The third advantage of the proposed form of the interaction is that the lattice equations can 
be used not only for the integer but also for fractional values of the parameter s. Therefore the 
suggested general lattice model can be extended on the fractional nonlocal case. The suggested 
types (1261) and (1281) of inter-particle interactions in the lattice can be used for non-integer s. 
If we consider interaction terms defined by (1261) and (1281) with non-integer s = a, then we 
will get continuum equations with the Riesz fractional derivatives [4T] of orders s = a by the 
methods suggested in PHESJ. The lattice models with long-range interactions of the types (1251) 
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and (1281) with non-integer s = a, can serve as microscopic models for elastic continuum with 
power-law non-locality. It allows us to derive fractional generalizations of gradient elasticity by 
using a microscopic approach [30, 3Tj. We also assume that the suggested lattice model can be 
generalized to get discrete (lattice) models for dislocations in the gradient elasticity continuum 
[321 [331, [31 [, (35[ 361 [37 , 381 09], and then it will be possible to extend them to the fractional 
non-local case. The suggested lattice models with long-range interactions, which are suggested 
for the gradient elasticity continuum, can be important to describe the non-local elasticity of 
materials at micro and nano scales [42, EES HI, S5j, where the interatomic and intermolecular 
interactions are prevalent in determining the properties of these materials. 


References 

[1] M. Born, K. Huang, Dynamical Theory of Crystal Lattices , (Oxford University Press, 
Oxford, 1954). 420 pages 

[2] A.A. Maradudin, E.W. Montroll, G.H. Weiss, Theory of Lattice Dynamics in the Harmonic 
Approximation (Academic Press, New York, 1963) 319 pages. 

[3] H. Botteger, Principles of the Theory of Lattice Dynamics (Academie-Verlag, Berlin, 1983) 
330 pages. 

[4] C. Kittel, Quantum Theory of Solids (Wiley, New York, London, 1963) 435 pages. 

[5] L.I. Sedov, A course in Continuum Mechanics Vol. I-IV. (Wolters-Noordhoff Publishing, 
Netherlands, 1971). 

[6] L.D. Landau, E.M. Lifshitz, Theory of Elasticity (Pergamon Press, Oxford, New York, 
1986) 187 pages. 

[7] E. Kroner, ” Elasticity theory of materials with long range cohesive forces”, International 
Journal of Solids and Structure. Vol.3. No.5. (1967) 731-742. 

[8] A.C. Eringen, ’’Nonlocal polar elastic continua”, International Journal of Engineering Sci¬ 
ence. Vol.10. No.l. (1972) 1-16. 

[9] A.C. Eringen, D.G.B. Edelen, ”On Nonlocal Elasticity”, International Journal of Engi¬ 
neering Science. Vol.10. No.3. (1972) 233-248. 

[10] A.C. Eringen, B.S. Kim, ’’Relation between non-local elasticity and lattice dynamics”, 
Crystal Lattice Defects. Vol.7. (1977) 51-57. 

[11] I.A. Kunin, Media with Microstructure Vol. I and II. (Springer-Verlag, Berlin, New York, 
1982 and 1983). 

[12] A.C. Eringen, Nonlocal Continuum Field Theories (Springer, New York, 2002) 376 pages. 

[13] D. Rogula, Nonlocal Theory of Material Media (Springer-Verlag, New York, 1983) 284 
pages. 

[14] H. Askes, E.C. Aifantis, ’’Gradient elasticity in statics and dynamics: An overview of 
formulations, length scale identification procedures, finite element implementations and 
new results”, International Journal of Solids and Structures. Vol.48. No.13. (2011) 1962- 
1990. 


15 


[15] H. Askes, A.S.J. Suiker, L.J. Sluys, ”A classification of higher-order strain-gradient models 
- linear analysis”, Archive of Applied Mechanics. Vol.72. No.2-3. (2002) 171-188. 

[16] R. Mindlin, ’'Theories of clastic continua and crystal lattice theories”, In: E. Kroner, 
(Eel.), IUTAM Symposium Mechanics of Generalized Continua (Springer-Verlag, Berlin, 
1968) pp. 312-320. 

[17] M. Rubin, P. Rosenau, O. Gottlieb, ”Continuum model of dispersion caused by an inherent 
material characteristic length”, Journal of Applied Physics. Vol.77. No.8. (1995) 4054-4063. 

[18] H.-B. Miihlhaus, F. Oka, ”Dispersion and wave propagation in discrete and continuous 
models for granular materials”, International Journal of Solids and Structures. Vol.33. 
No.19. (1996) 2841-2858. 

[19] A.V. Metrikine, H. Askes, ” One-dimensional dynamically consistent gradient elasticity 
models derived from a discrete microstructure. Part 1: Generic formulation”, European 
Journal of Mechanics - A/Solids. Vol.21. No.4. (2002) 555-572. 

[20] W. Chen, J. Fish, ”A dispersive model for wave propagation in periodic heterogeneous 
media based on homogenization with multiple spatial and temporal scales”, Journal of 
Applied Mechanics. Transactions ASME. Vo.68. No.2. (2001) 153-161. 

[21] H. Askes, A. Metrikine, ”One-dimensional dynamically consistent gradient elasticity mod¬ 
els derived from a discrete microstructure. Part 2: Static and dynamic response”, European 
Journal of Mechanics. A. Solids. Vol.21. No.4. (2002) 573-588. 

[22] H. Askes, A. Metrikine, ’’Higher-order continua derived from discrete media: continuali- 
sation aspects and boundary conditions”, International Journal of Solids and Structures. 
Vol.42. No.l. (2005) 187-202. 

[23] V.P. Maslov, Operator Methods (Mir, Moscow, 1976) or (Nauka, Moscow, 1973 in Russian), 
Section 8. 

[24] V.E. Tarasov, ’’Map of discrete system into continuous”, Journal of Mathematical Physics. 
Vol.47. No.9. (2006) 092901. (arXiv:0711.2612) 

[25] V.E. Tarasov, ’’Continuous limit of discrete systems with long-range interaction”, Journal 
of Physics A. Vol.39. No.48. (2006) 14895-14910. (arXiv:0711.0826) 

[26] V.E. Tarasov, G.M. Zaslavsky, ’’Fractional dynamics of coupled oscillators with long-range 
interaction”, Chaos. Vol.16. No.2. (2006) 023110. (nlin.PS/0512013) 

[27] V.E. Tarasov, Fractional Dynamics: Applications of Fractional Calculus to Dynamics of 
Particles, Fields and Media (Springer, New York, 2011) 504 pages. 

[28] A.P. Prudnikov, Yu.A. Brychkov, O.I. Marichev, Integrals and Series, Vol. 1: Elementary 
Functions (Gordon and Breach, New York, 1986) 

[29] H. Askes, L.J. Sluys, ”A classification of higher-order strain-gradient models in damage 
mechanics”, Archive of Applied Mechanics. Vol.73. No.5. (2003) 448-465. 

[30] V.E. Tarasov, ’’Lattice model with power-law spatial dispersion for fractional elasticity”, 
Central European Journal of Physics. Vol.11. No.11. (2013) 1580-1588. (arXiv:1501.01201) 

[31] V.E. Tarasov, ’’Lattice model of fractional gradient and integral elasticity: Long-range 
interaction of Grunwald-Letnikov-Riesz type”, Mechanics of Materials. Vol.70. No.l. (2014) 
106-114. 


16 



[32] M.Yu. Gutkin, E.C. Aifantis, ’’Dislocations and disclinations in the gradient theory of 
elasticity”, Physics of the Solid State. Vol.41. No.12. (1999) 1980-1988. 

[33] M.Yu. Gutkin, E.C. Aifantis, ’’Dislocations and disclinations in gradient elasticity”, Phys- 
ica Status Solidi B. Vol.214. No.2. (1999) 245-284. 

[34] M. Lazar, G.A. Maugin, ’’Nonsingular stress and strain fields of dislocations and discli¬ 
nations in first strain gradient elasticity”, International Journal of Engineering Science. 
Vol.43. No.13-14. (2005) 1157-1184. (arXiv:cond-mat/0502023) 

[35] M. Lazar, G.A. Maugin, E.C. Aifantis, ”On dislocations in a special class of general¬ 
ized elasticity”, Physica Status Solidi B. Vol.242. No.12. (2005) 2365-2390. (arXiv:cond- 
mat/0504291) 

[36] M. Lazar, G.A. Maugin, E.C. Aifantis, ’’Dislocations in second strain gradient elasticity”, 
International Journal of Solids and Structures. Vol.43. No.6. (2006) 1787-1817. 

[37] M. Lazar, G.A. Maugin, ’’Dislocations in gradient elasticity revisited”, Proceedings of the 
Royal Society A. Vol.462. No.2075. (2006) 3465-3480. 

[38] M. Lazar, ’’Non-singular dislocation loops in gradient elasticity”, Physics Letters A. 
Vol.376. No.21. (2012) 1757-1758. (arXiv: 1204.0945) 

[39] M. Lazar, ’’The fundamentals of non-singular dislocations in the theory of gradient elas¬ 
ticity: Dislocation loops and straight dislocations”, International Journal of Solids and 
Structures. Vol.50. No.2. (2013) 352-362. (arXiv: 1209.1997) 

[40] A. Erdelyi, W. Magnus, F. Oberhettinger, F.G. Tricomi, Higher Transcendental Functions 
Volumes 1. (McGraw-Hill, New York, 1953), or (Krieeger, Melbourne, Florida, 1981). 

[41] A.A. Kilbas, H.M. Srivastava, J. J. Trujillo, Theory and Applications of Fractional Differ¬ 
ential Equations (Elsevier, Amsterdam, 2006) 540 pages. 

[42] M. Ostoja-Starzewski, ’’Lattice models in micromechanics”, Applied Mechanics Reviews. 
Vol.55. No.l. (2002) 35-60. 

[43] L. Merchan, R. Szoszkiewicz, E. Riedo, ” Nano Mechanics: elasticity in nano-objects”, in 
Fundamentals of Friction and Wear NanoScience and Technology. E. Gnecco, E. Meyer 
(Eds.), (Springer, 2007) pp.219-254. 

[44] E.C. Aifantis, ’’Gradient effects at macro, micro, and nano scales”, Journal of the Mechan¬ 
ical Behavior of Materials. Vol.5. No.3. (1994) 355-375. 

[45] E.C. Aifantis, ’’Gradient nanomechanics: applications to deformation, fracture, and dif¬ 
fusion in nanopolycrystals”, Metallurgical and Materials Transactions A. Vol.42. No. 10. 
(2011) 2985-2998. 


17 




Figure 2: Plot of the function (l27j) for the range x G [0,4] and y G [0,6]. 



Figure 3: Plot of the function (1271) for the range x G [2,4] and y G [1,3]. 
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